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ABSTRACT 

Quasi-constant heating at the footpoints of loops leads to evaporation and condensation cycles of the 
plasma; thermal non-equilibrium (TNE). This phenomenon is believed to play a role in the formation of 
prominences and coronal rain. However, it is often discarded to be involved in the heating of warm loops 
as the models do not reproduce observations. Recent simulations have shown that these inconsistencies 
with observations may be due to oversimplifications of the geometries of the models. In addition, our 
recent observations reveal that long-period intensity pulsations (several hours) are common in solar coro¬ 
nal loops. These periods are consistent with those expected from TNE. The aim of this paper is to derive 
characteristic physical properties of the plasma for some of these events to test the potential role of TNE 
in loop heating. We analyzed three events in detail using the six EUV coronal channels of SDO/AIA. We 
performed both a Differential Emission Measure (DEM) and a time-lag analysis, including a new method 
to isolate the relevant signal from the foreground and background emission. Eor the three events, the 
DEM undergoes long-period pulsations, which is a signature of periodic heating even though the loops 
are captured in their cooling phase, as is the bulk of the active regions. We link long-period intensity pul¬ 
sations to new signatures of loop heating with strong evidence for evaporation and condensation cycles. 
We thus witness simultaneously widespread cooling and TNE. Einally, we discuss the implications of our 
new observations for both static and impulsive heating models. 

Subject headings: Sun: corona - Sun: oscillations - Sun: UV radiation 


1. Introduction 


The heating mechanism(s) of coronal loops that 
generate million-degree plasma and maintain them at 
such temperatures remain unknown. Some of the 
main challenges is to match observational results with 
model predictions and to find signatures specific to 
some heating process. Among heating models, sce¬ 
narios of two different kinds can be differentiated: the 
impulsive heating scenarios and the quasi-continuous 
ones. 


Small-scale impulsive heating scenarios are known 
under the term of ’’nanoflares”. A nanoflare consists 
in a relatively small energy impulsive event occurring 


by reconnection of loop strands (e.g Cargill & Klim- 


|chuk|199Tl|Patsourakos & Klimchuk|20061 |Klimchuk| 


20061. If many of these stochastic events happen close 
in time in the loops, we have what is called a nanoflare 
storm. The nanoflare term was first proposed by|Parker 


( |1988| l, and, since then, signatures of such events are 
investigated. The heating by nanoflares can produce 
million-degree type loops ( |Guarrasi et al.||2010[ ) as 
well as very high temperature loops {T >5 MK) 
( Cargill|T994[ Klimehuk 20061, which is one of the 
main strong points of heating by nanoflares. 


The other class of scenarios involve continuous 
or quasi-continuous heating concentrated at the foot- 
points of the loops. Numerical simulations show that 
in that case, the loops can undergo cycles of evapora¬ 
tion and condensation, a phenomenon called thermal 
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non-equilibirum. Thermal non-equilibirum processes 
are known to play an important role for prominences 


(Antiochos & Klimchuk|19^ Karpen et al. 20061 and 


coronal rain ( [Muller et al.|2005 1. Coronal rain is con¬ 
sidered to be the consequence of ’’catastrophic cool¬ 
ing” events, as seen in numerical simulations of ther¬ 
mal non-equilibrium ( |Muller et al.||2003l |2004| ). |An-| 
[tolin et aklpOlOl l connect coronal rain to coronal heat¬ 
ing mechanisms by means of a comparison between 
observations and MHD simulations of loops heated at 
the footpoints. This heating produces condensation- 
evaporation cycles that, in case of a strong plasma 
condensation at the loop apex, form ’’blobs” of cold 
plasma {T ~ 0.1 MK). Motion of these blobs were 
observed by De Groof et al. ( 2004|l with the Extreme 


Ultraviolet Imaging Telescope (EIT; Delaboudiniere 


et al.||1995|l on board the Solar and Heliospheric Ob¬ 


servatory (SoHO; [Domingo et al. 11993) 1 and jSchrijverj 
([200 l|l with the Transition Region And Coronal Ex¬ 
plorer (TRACE; |Handy et al.|1999l l. 

|Mok et al.| P008| l consider that thermal non¬ 
equilibrium could be involved also in the formation 
of coronal loops. In their simulations, a heating 
concentrated at the footpoints leads to condensation- 
evaporation cycles. The temperature oscillates in time 
with a period of about 16 hours. However, the role 
of thermal non-equilibrium in heating of warm coro¬ 
nal loops has often been discarded, as simulations are 
not able to reproduce all the observational constrains 
( jKlimchuk et al.|2010| l. 

Nevertheless, oversimplifications of the model ge¬ 
ometry could explain these inconsistencies with obser¬ 
vations. jLionello et al.| ( |2OT3| l show an analysis of syn¬ 
thetic images of an active region, obtained from 3D- 
hydrodynamic simulations, that seems consistent with 
observational results. In their paper, the authors ana¬ 
lyzed emission from loops in the same way as if they 
analyzed real images with diagnostics confirmed by 


Winebarger et al. (20141. They successfully checked 


that the loops extracted reveal seven fundamental char¬ 
acteristics of warm loops. In this series of papers. 


Mikic et al. (2013 i explore influences of the geom¬ 


etry of the loops and of the symmetry of heating in 
ID models. Some cases lead to ’’catastrophic cool¬ 
ing” events that they called ’’complete” condensations. 
However, with other configurations loops stay stable 
to temperatures above 1 MK when only ’’incomplete” 
condensations are formed. 


Recently, long-period (several hours) intensity pul¬ 
sations have been found to be very common in ac¬ 
tive regions, especially in loops ( Auchere et al.|20f^ . 
These pulsations are likely to be new observational sig¬ 
natures of the heating processes in coronal loops. 

In this paper we analyze in detail some long-period 
pulsation events found in loops to investigate the phys¬ 
ical mechanisms involved, using data from the Atmo¬ 
spheric Imaging Assembly (AIA; Boerner et al.|2012 


[Lemen et al.|2012^ on board the Solar Dynamics Ob¬ 
servatory (SDO; Pesnell et al.||2012] l. This allows to 
track long-period intensity pulsation events simultane¬ 
ously in six coronal extreme ultraviolet (EUV) chan¬ 
nels; 94 A, 131 A, 171 A, 193 A, 211 A, 335 A. 
Compared to the previous analysis in one EIT band 
( Auchere et al.|20f^ , we can thus now perform phys¬ 
ical diagnostics of the plasma. The three cases that 
we present here show strong evidence for cycles of 
evaporation and of incomplete condensation. Eirst we 
present the analyzed data in Sec. 2.1 Then we move 


to the physical analysis with a DEM diagnostic in 


Sec. 2.2 and to giving evidence for widespread cool¬ 
ing in Sec.|2.3| 


2. Analysis of 3 typical loop events 
2.1. Data sample 

In order to detect long-period intensity pulsation 
events, we use the automatic detection algorithm de¬ 


veloped by Auchere et al. (20141. The algorithm treats 
sequences of about 6 days of data at a cadence of 13 
minutes, which represents about 600 images. Eor each 
sequence, a region of interest (ROI) is tracked by map¬ 
ping these images into heliographic coordinates (with 
a resolution of 0.2° in longitude and latitude) and by 
compensating for the differential rotation. The data 
cube is analyzed in the Eourier space (obtained by 
Eourier transform of the data cube along the time axis) 
with respect to several criteria, the main one being a 
detection threshold at lOcr above an estimate of the 
average local power. Using the database of images of 
EIT/SoHO in the 195 A passband, from January 1997 
to July 2010, Auchere et al. ( |2014| ) reported 917 events 
of long-period intensity pulsations. About half of these 
events have been detected in active regions and 50% of 
these have been visually associated with loops. 

We analyze here data from AIA using the same de¬ 
tection algorithm. The AIA raw data are read with 
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Fig. 1.— Individual frames for event 1 for the six coronal channels of SDO/AIA on June 07, 2012 at 13:47 UT. 
This event is localized in NOAA AR 11499. These images are displayed with a logarithmic scale and in heliographic 
coordinates. The green contour delimits the area of the pulsations detected at 335 A. The black contour delimits the 
area manually selected for detailed time series analysis. 
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the routine read_sdo from the Interactive Data Lan¬ 
guage SolarSoftware library. We normalize the inten¬ 
sities by the exposure time and so they are expressed 
in DN.pix *.s ' (DN: Digital Number). In order to in¬ 
crease the signal-to-noise ratio, we start from 4x4 
binned images since the detected regions are large 
(larger than the typical area of bright points). In ad¬ 
dition, given the detected periods, we use a cadence 
of 13 min. The differential rotation is compensated 
for all the passbands wi th the rate measured in the 
195 A passband of EH ^Hortin||2003 Auchere et al. 


|2014[ ), which corresponds to the average 1 MK corona. 
For more details about the detection algorithm, the 
readers should refer to [Auchere et al. P014|l. For 


the Fourier analysis (sections 2.1 and|2.2|) and all the 


cross-correlation computations, we resample the time 
series using a linear interpolation to ensure a regular 
cadence. We checked that the resampling does not 
affect the detections of long-period pulsations in the 
Fourier analysis, as it affects only high frequencies. 
See |Auchere et al. (2014 Appendix A: Sources of spu¬ 
rious frequencies) for details about sources of spuri¬ 
ous detections. Since the time series have about 300 
images (sequences of three days of data), the lOcr de¬ 
tection threshold corresponds to a confidence level of 
99%. 


From May 2010 to December 2013 the code de¬ 
tected about 2000 events. As for the FIT based study, 
50% of them are localized in active regions and about 
a half of them are visually associated with loops. Since 
we detect some events that have also been detected 
with FIT (in May 2010), we can discard definitely in¬ 
strumental artefacts as the cause for these long-period 
intensity pulsations. Furthermore, since we do not ap¬ 
ply any filter to the data before Fourier analysis, it 
is also unlikely that these events are an artefact of 
processing. Among the detected events, we choose 
to analyze three events showing a strong signal in 
loops, in order to investigate the underlying physi¬ 
cal processes in detail. These three events cover a 
large range of periods. For the first event, the power 
peaks around 30.7 //Hz, i.e. 9.0 hours, for the sec¬ 
ond around 49.9 //Hz, i.e. 5.6 hours and for the third 
around 72.9 //Hz, i.e. 3.8 hours. We can thus explore 
the possible relationship between frequency and phys¬ 
ical properties. The three active regions with pulsating 


loops are also very different with respect to their size. 
In addition, the three tracked active regions have a 
moderate activity with no flares above C-class. To sim¬ 
plify the reading and because we reach the same con¬ 
clusions with the three events, we only present event 1 
in detail in the core of the text. Events 2 and 3 are 
presented respectively in Appendix |A| and AppendixjB] 
with the same progression. 


Fig. represents the heliographic held of view for 
event 1 on June 07, 2012 at 13:47 UT in the six coronal 
channels of SDO/AIA. We tracked NOAA AR 1499 
during 154 hours (i.e. more than 6 days), from June 
03, 2012 18:00 UT to June 10, 2012 04:29 UT. In order 
to minimize distortion effects due to the heliographic 
mapping transformatiorj^ we kept only three central 
days in this sequence of six days, from June 05, 2012 
11:14 UT to June 08, 2012 11:16 UT. For the rest of the 
data analysis we will use this short sequence, except if 
otherwise stated. The green contour in Fig. [T] delimits 
the area of pulsations detected at 335 A, which is the 
passband with the strongest signal (see. Fig. |^. This 
contour is centered on 182.9° of longitude and 25.3° 
of latitude and has an area of 45.5 heliographic square 
degrees (corresponding to 5271 Mm^ on the sphere). 
This contour fits the shape of large loops visible in 
images. We manually selected a smaller contour in 
black to delimit the area used for detailed time series 
analysis. This contour is included in the contour de¬ 
tected at 335 A and the signal is strong for most of the 
passbands: more than Icr of normalized power, aver¬ 
aged over this contour for five passbands (see Fig. [^. 
We choose this contour so that the phases of the pe¬ 
riodic components (see Fig. 10 1 are similar in order 
to maximize the signal-to-noise ratio when averaging 
over severals heliographic pixels. 


We choose to concentrate on time series of the av¬ 
erage parameters in this black contour, because isolat¬ 
ing pulsating loops and subtracting a background (al¬ 
beit crucial for most studies of coronal loops) is nei¬ 
ther practical nor meaningful for the present analy¬ 
sis. Actually, our time series are much longer than the 


typical life time of loops: 1000-5000 s ( 

Winebarger 

et al. 

2003 

Winebarger & Warren||2005[ 

Jgarte-Urra 

et al. 

2009 

1 , therefore it is not possible to isolate and 


to track loops over such long times. In Fig.|^ we dis¬ 
play the six light curves averaged over this black con- 


' We assume that the structures rotate with the same velocity indepen- 
dently of the channel used. 


^Distortion of the stnactures is more important close to the limb and 
if the stnactures ai‘e tracked for a long time. 


4 
































60 

50 

40 

30 

20 

10 

0 


emission measure — peak temperature — slope — 



j_^_I_I_I_I_I_U_I_I_I_I_I_I_I_I_I_I_I_I_I_I_ Li_I_I_I_I_I_^_L 

20 40 60 80 100 120 140 

Date (hours since 2012-06-03 18:00) 


5 




Fourier power spectrum (94 A) 



Fourier power spectrum (171 A) 



Fourier power spectrum (211 A) 



Frequency (qHz) 


Fourier power spectrum (131 A) 



Fourier power spectrum (193 A) 



Fourier power spectrum (335 A) 
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Fig. 3.— Fourier power spectra for the six AIA coronal channels. These spectra are computed for the light curves 
between the two vertical dashed lines in Fig.[^ The dashed line is the estimate of the average local power. The solid 
line is the lOcr detection threshold. If a peak of power is above the detection threshold, we display the corresponding 
period of pulsations and normalized power at the top right of each spectrum. 
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toured area: 94 A in black, 131 A in blue, 171 A in 
cyan, 193 A in green, 211 A in orange and 335 A in 
red. These light curves are normalized to standard de¬ 
viation (we subtract the mean curve and divide by the 
standard deviation). The vertical dashed lines delimit 
the short sequence (i.e. 3 days) used to compute the 
Fourier spectra. 


Fig. [^represents the corresponding Fourier power 
spectra, normalized to the variance cTq of these light 
curves ( Torrence & Com^|1998| l. The dashed line 
is the estimate of the average local power. The solid 
line is the lOcr detection thresholcQ We consider only 
peaks of power between 18 pHz and 110 fiHz (be¬ 
tween 15.4 and 2.5 hours). In the 94 A, 131 A, 193 A, 
211 A, and 335 A channels there is a significant peak 
of power at 30.7 fjHz (9.0 hours). The normalized 
power at these peaks are respectively 12.9cr, 10.6cr, 
17.5cr, 22.4cr, and 40.9cr. In the 171 A passband there 
is no peak of normalized power above lOcr. However, 
there is a 8.9cr peak at the frequency detected in the 
other bands (30.7 fjHz), plus a second 9.3cr peak at 
46.1 fiHz (6.0 hours). This small peak at 6.0 hours ex¬ 
ists also in the 94, 131, 211, and the 335 channels with 
respectively 7.4cr, 8.3cr, 8.0cr, and 8.9cr of normalized 
power. All these peaks correspond to confidence lev¬ 
els greater than 90%. At 211 A and 335 A, we can also 
notice small peaks at 42.2 //Hz with respectively 7.4cr 
and 9.2cr of normalized power. 


Fig. [^represents the maps of normalized power in 
heliographic coordinates for the six passbands at the 
frequency of detection (30.7 //Hz, i.e. 9.0 hours), with 
the same field of view in Fig.[2 The normalized power 
is displayed in logarithmic scale, with a saturation at 
20cr. As for Fig. [T] we displayed the contour detected 
at 335 A in green and the contour that we use for anal¬ 
ysis in black. We can notice that the power is well lo¬ 
calized in loops for most of the passbands. In the black 
contour, the normalized power is higher or equal to lo¬ 
on average for all the passbands, except for 131 A the 
power is the weakest (5cr in average). This value rises 
to 19cr at 335 A. 


^cr is the averaged power as a function of frequency and is different 
from cTg, which is the variance of the time series. In the case of 
white noise, we would have cr = cr^ 


2.2. Differential emission measure (DEM) analy¬ 
sis 

2.2.1. Method 


We compute the differential emission measure 


(DEM) using the method developed in Guennou et al. 


( 2012a|b| l. To fit the DEM, we choose the active region 
model developed by Guennou et al. (2013 i with the 
six coronal EUV passbands of AIA. This parametric 
model is designed to represent the shape of the DEM 


usually observed in active regions (Warren et al. 2011 


Winebarger et al. 201 l| l. Eour parameters determine 
the shape (as in Eig.[7]l of this DEM model: 


• the slope a of the low temperature wing repre¬ 
sented by a power law 

• the peak temperature of the DEM 

• the total emission measure (the integral of the 
DEM over the electron temperature T^) 

• the width cr of the high temperature Gaussian 
wing. We choose to set this parameter to a fixed 
value of 0.1 log T as the high temperature wing 
of the DEM is poorly constrained by AIA data. 


We compute the DEM using the six channels of 
AIA. Eor this, we grouped images from the six pass- 
bands into sextuplets allowing for a maximum non si¬ 
multaneity of +2 minutes. 


2.2.2. Results 

Eor event 1, the evolution of the slope, peak tem¬ 
perature, emission measure, and residuals (a mea¬ 
sure of the adequacy of the chosen model to represent 
the real DEM), averaged over the black contour intro¬ 
duced in Eig.[2 is presented in Eig.[^ Eig. [^represents 
the power spectra of these time series. The peak tem¬ 
perature, the emission measure and x^ present pulsa¬ 
tions at the same frequency as the intensity. Actually 
there are power peaks at 30.7 //Hz (9.0 hours) with re¬ 
spectively 20.1cr, 19.9cr and 12.8cr for these two DEM 
model parameters and for x^. Eor the slope there is 
no significant power peak but a small peak (7.4cr) at 
46.1 //Hz, the frequency of the secondary peaks found 
in the power spectra of the light curves (see Eig. [^. 
Even if we found pulsations in the absolute val¬ 
ues of x^ varie between 0.2 and 1.5, which indicates a 
generally good fit to the DEM. 
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Fig. 4.— Maps of normalized power at 30.7 pHz (i.e. 9.0 hours) for event 1. These maps are in heliographic 
coordinates for the six EUV passbands of AIA: 94 A, 131 A, 171 A, 193 A, 211 A, 335 A. The normalized power is 
in logarithmic scale, from 1 cr to 20cr. The black contour, included in the detected contour in green, is the one used for 
detailed time series analysis. 
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Fourier power spectrum (siope) 
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Fourier power spectrum (emission measure) 



Fourier power spectrum (x^) 
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Fig. 5.— Fourier power spectra for the three parameters of the DEM model for event 1: the peak temperature, the 
emission measure and the slope. We represent also the Fourier power spectrum for^^. These spectra are computed 
for light curves between the two vertical dashed lines in Fig. The dashed line is the estimate of the average local 
power. The solid line is the lOcr detection threshold. If a peak of power is above the detection threshold, we display 
the corresponding period of pulsations and normalized power at the top right of each spectrum. 
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Time shift (minutes) 

Fig. 6.— Cross-correlation values for peak tempera¬ 
ture/slope in solid line and peak temperature/emission 
measure in dotted line. We explored time shifts from 
-300 minutes to 300 minutes. Time lag for each pair of 
DEM parameters is indicated by a black dot. We used 
the averaged curves at 1 min of cadence averaged over 
the small black contour presented in Fig. [T] between 
June 06, 2012 08:42 UT (63 h after the beginning of 
the long sequence) and June 07, 2012 12:00 UT (90 h 
after the beginning of the long sequence). 

Although there is no significant power at the fre¬ 
quency detected in the intensity time series for the 
DEM slope, we found a correlation between the slope 
and both the emission measure and the temperature, as 
seen in Fig. The cross-correlation between the peak 
temperature and the emission measure is in dotted line 
and that between the peak temperature and the slope is 
in solid line with a time shift between -300 min and 
-1-300 min. These curves are computed over a duration 
of about 28 hours, between 62.7 hours after the begin¬ 
ning of the sequence (i.e. on June 06, 2012 08:42 UT) 
and 90 hours after the beginning of the sequence (i.e. 
on June 07, 2012 12:00 UT). There is a time lag of -40 
minutes for a peak cross-correlation value of 0.7. This 
means that variations of the slope precede variations 
of the peak temperature. Between the peak tempera¬ 
ture and the emission measure, there is also a correla¬ 
tion. The peak temperature is in advance of the emis¬ 
sion measure with a time lag of 119 minutes for a peak 
cross-correlation value of 0.5. 

If we look at the variations of the absolute values of 
the DEM parameters (instead of normalizing them by 
their standard deviation, as in Fig. |^, we can notice 



Fig. 7.— Four extreme cases of the shape of the DEM 
model (see Fig. |^. In green, the shape of the DEM 
on June 05, 2012 11:14 UT with the minimum slope 
(case 1) and in red on June 06, 2012 21:42 UT with the 
maximum slope (case 3). In orange, the shape of the 
DEM on June 06, 2012 08:42 UT with the minimum 
peak temperature (case 2) and in blue on June 07, 2012 
15:30 UT with the maximum peak temperature (case 
4). The four dates corresponding to these cases are 
represented by plus signs in Fig. 

that the amplitude of the variations of the peak tem¬ 
perature is relatively small compared to the variations 
of the emission measure and of the slope: the peak 
temperature varies from 2.8 MK to 3.1 MK (a rela¬ 
tive amplitude of 11%), while there are variations from 
2.03 X 10^^ cm“^ to 2.82 x 10^’ cm“^ for the emission 
measure (relative amplitude of 38%) and from 2.7 to 
4.7 for the slope (relative amplitude of 71%). We can 
note that larger temperature variations are found for 
event 2 and event 3 (see Appendix [A| and [B]l . 

Fig. represents four extreme cases of the shape 
of the DEM during the sequence. Instants chosen for 
the four cases are pointed out sequentially by black 
plus signs in Fig. This figure illustrates that, of 
the three parameters of the DEM model, the slope 
undergoes the largest relative variations. In case 1 
(in green), the DEM model is the most multither- 
mal, as the slope reaches a minimum, a - 2.7 (on 
June 05, 2012 11:14 UT), the total emission measure 
is 2.49 X 10^^ cm“^ and the peak temperature is at 
3.0 MK. In orange (case 2), we trace the shape of 
the DEM model when the peak temperature reaches 
a minimum at 2.8 MK. That happens on June 06, 


10 




2012 08:42 UT when the slope a - 4.0 and the total 
emission measure is 2.29 x 10^^ cm“^. In red (case 3), 
on June 06, 2012 21:42 UT, the slope is the steepest 
with a = 4.7. The peak temperature is at 2.9 MK and 
the total emission measure is 2.67 x 10^^ cm“^. For the 
case 4 (in blue), the peak temperature reaches a maxi¬ 
mum at 3.1 MK. The slope is weaker with a = 3.0 and 
the total emission measure is 2.82 x 10^^ cm“^. This 
last case occurs on June 07, 2012 15:30 UT. 


These results show that the DEM changes from a 
multithermal distribution to a more isothermal distri¬ 
bution and that at least the peak temperature and emis¬ 
sion measure vary periodically. In addition, even if 
the normalized power found in the DEM slope Eourier 
spectrum is under the detection threshold, there is a 
clear correlation between slope and peak temperature, 
which is another indication that the slope also under¬ 
goes periodic variations. Therefore, we conclude from 
the DEM and correlation analysis that there are peri¬ 
odical changes in the thermal structure of the loops 
where we have detected long-period intensity pulsa¬ 
tions. These periodical modifications of the ther¬ 
mal structure must be then a consequence of heat¬ 
ing in loops, given that a cooling phase is necessar¬ 
ily preceded by a heating phase. This indicates that 
these long-period intensity pulsations are linked to 
loop heating. 


2.3. Evidence for widespread cooling 

2.3.1. Method 


We have seen that the thermal structure of the pul¬ 
sating loops undergo periodical changes which links 
these phenomena to loop heating. However when 
the thermal structure of these loops evolves from a 
multithermal case to a more isothermal case, we do 
not observe heating. As was shown by several au¬ 
thors, EUV loops are generally seen in their cool¬ 
ing phase ( Winebarger&^Wanenjj2005 Warren et al. 


2002 [Winebarger et al.|2003[ Ugarte-Uira et al.|2006 


2009 [Mulu-Moore et aL|2011| l' In order to show 

that this is the case for a large fraction of the active 
region in which pulsations are detected, we use the 
method developed in Viall & KIimchuk| (20I2 here¬ 
after V&K12). V&K12 present an analysis of an ac¬ 
tive region and show that the plasma is mainly in state 
of cooling, as the EUV intensity peaks first in the hot¬ 
ter passbands and then in the cooler passbands. By 
computing the time lag between several pairs of chan¬ 


nels, V&K12 obtain the temporal succession of chan¬ 
nels and see that the plasma is cooling down. By using 
2 hour and 12 hour time series, V&K12 derived time 
lags ranging from a couple of minutes to 1.5 hours. 
Since the lines dominating the passbands are formed 
by collisions, variations of the emission measure pro¬ 
duce simultaneous variations in the six coronal pass- 
bands. On the contrary temperature variations are ex¬ 
pected to be gradually reflected in the passbands ac¬ 
cording to the ordering of peak temperature response 
of the channels: 335 A (2.5 MK), 211 A (2 MK), 
193 A (1.5 MK), 94 A (1 MK and 7 MK), 171 A (0.8 
MK) and 131 A (0.5 MK) (AIA temperature response 
functions; [Boerner et al.|2012{ |Lemen et al.|2012] l. 

2.3.2. Results 



Time shift (minutes) 


Eig. 8.— Cross-correlation values for six pairs of AIA 
channels (event 1): 94-335 (red), 335-211 (blue), 211- 
193 (green), 335-193 (orange), 335-171 (cyan) and 
171-131(black). We used light curves at 1 min of ca¬ 
dence averaged over the small black contour presented 
in Eig. [T] We explored time shifts from -300 minutes 
to 300 minutes. The time lag for each pair of channel 
is indicated by a colored dot. 


Our results in Eig. are to be compared with Eig. 4 
of V&K12. As detailed below, we observe the same 
behavior as them, meaning that the plasma in the small 
black contoured area is cooling. 

Eig. [^represents the cross-correlation values for six 
pairs of AIA channels: 94-335, 335-211, 211-193, 
335-193, 335-171, and 171-131. This calculation is 
made with the averaged light curves in the small black 
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contour seen in Fig.j^ Here we choose a time cadence 
of 1 min, as we expect time lags under the 13 min ca¬ 
dence used previously. We explore time shifts from 
-300 minutes to 300 minutes (i.e. 5 hours). Time lags 
are given by the peak cross-correlation value for each 
pair of channels. 

For the 94-335 pair, we obtain a time lag of -121 
minutes. This negative time lag indicates that the in¬ 
tensity peaks in the 335 channel before the 94 chan¬ 
nel. The temperature response of the 94 channel has 
two peaks, one cooler than the 335 peak temperature 
and one hotter. Since 94 peaks after 335, we conclude 
that the plasma does not reach the temperature of the 
hot peak of the 94 band. The 193 channel follows the 
211 channel with a time lag of 30 min. For the pairs 
335-211, 335-193, and 335-171, time lags are positive 
and increase gradually, with 115 min, 144 min and 153 
min respectively. The 171 and 131 channels show no 
significant time lag (-1 min) given the time resolution. 

Conclusions are more difficult to draw for the 94- 
335, 335-193, and 335-171 pairs, which peak cross¬ 
correlation values are respectively 0.3, 0.4, and 0.3. 
Peaks for 335-193 and 335-171 cross-correlations are 
not well separated and the time lag for 94-335 is 
smaller in absolute value than the time lag between 
335-193 and 335-171. It could be explained because 
the 94, 171 and 193 channels are the most distant in 
time compared to the 335 channel; the emission mea¬ 
sure is more likely to change during a longer time and 
the intensities in these passbands are thus more likely 
to be de-correlated from that in 335. Moreover, we use 
longer time series that in V&K12 and thus the pass- 
bands are more likely to be de-correlated. Superimpo¬ 
sition of structures along the line of sight could also 
explain this effect since the background contribution 
may be different in the different passbands. 

This analysis can be generalized to the whole ac¬ 
tive region. We build time lag maps like Fig. 5 to 7 
of V&K12. For each pair of channels, the time lag for 
each pixel is given by the peak cross-correlation value 
of the intensity time series. Our time lag analysis is 
fully consistent with the widespread cooling observed 
by V&K12. In particular, we observe the same order¬ 
ing of channels. We can notice that this behavior is 
found in pulsating areas like in the rest of the active 
region. 

Fig. 1^ represents time lag maps for the entire ac¬ 
tive region. Black areas mark regions where the peak 


cross-correlation value is under 0.2. In red, the time 
lag is positive, in blue negative and in white we have 
a zero time lag. The 335-211, 211-193, 335-193, 335- 
171 maps are dominated by positive time lags. In the 
171-131 map we mainly observe a zero time lag. The 
94-335 map reveals two different behaviors. In the 
core of the active region, the time lag is mainly pos¬ 
itive, which indicates that the high temperature com¬ 
ponent of 94 dominates. On the contrary in the outer 
loops (including the black contour) the time lag is 
mainly negative: the low temperature component of 
94 dominates. With this method, we do not try to sub¬ 
tract the background and the foreground to isolate the 
loops. Thus, superimpositions of different structures 
along the line of sight can occur. 

To conclude on the behavior of the pulsating loops 
with respect to the widespread cooling of the active 
region, we isolate the pulsating structure in Fourier 
space; we subtract phase maps at the detected pulsa¬ 
tion frequency for pairs of channels, and we convert 
the phase difference into a time lag. Fig.[T^shows the 
phase maps. There is no clear evidence for variation 
of phases along the loops. Fig. [^displays the time 
lag maps obtained with the method of phase maps dif¬ 
ference. As can be seen in these maps, the phase dif¬ 
ference is only meaningful in areas with detected pul¬ 
sations; in other regions, the phases are random and 
so are the time lags. With this method, the pulsat¬ 
ing loops have been isolated, and the background and 
foreground emission, which are assumed not to be pul¬ 
sating at the detected frequency, are not taken into ac¬ 
count in the calculation of time lags. In addition, we 
can discern the same pattern of time lags in the pul¬ 
sating loops as in Fig. which was obtained with the 
peak cross-correlation method. To use phase maps to 
recover time lag between channels allows to discard 
contribution of background and foreground. More¬ 
over, we extract the signal of the pulsating loops pixel 
by pixel without any assumption on the shape of the 
structure. 

We compare time lags averaged over the contour ob¬ 
tained with the peak cross-correlation method and with 
the phase difference method in Table. With both 
methods we find similar time lags and approximately 
the same ordering of channels. For method 1 (cross¬ 
correlations), the time lag for 335-193 is larger than 
the time lag for 335-171 (145 min compared to 142 
min). So the intensity peaks first in the 171 channel 
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Fig. 9.— Time lag maps made with peak cross-correlation values for the same six pairs of AIA channels as in Fig.|^ 
94-335, 335-211, 211-193, 335-193, 335-171 and 171-131. We explored time shifts from -300 minutes to 300 minutes. 
Black areas represent peak cross-correlation values under 0.2. The green contour is the detected contour and the black 
contour is the contour chosen to trace the light curves. 


Table 1; Comparison of time lag values between the two methods: peak cross-correlation values and differences of 
phase. These are the averaged time lag in minutes in the black contour. 


Pairs of channels Time lag (min) Time lag (min) 

from cross-correlations (Fig.[^ from differences of phase Fig .(11 1 


335-211 

113 

122 

211-193 

26 

20 

335-193 

145 

142 

94-335 

-115 

-114 

335-171 

142 

155 

171-131 

-1 

-51 
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Fig. 10.— Phase maps at the frequency of detection (30.7 yuHz, i.e. 9.0 hours) for event 1. The phase is displayed in 
minutes with a cyclic color table. Here the detected contour is in black to facilitate the reading. Areas with normalized 
power under 4cr are in white. 
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Fig. 11.— Time lag maps obtained from differences of phase maps (see Fig.[T0)l for the same six pairs of AIA channels 
as in Fig.|g 94-335, 335-211, 211-193, 335-193, 335-171, and 171-131. We explored time shifts from -300 minutes 
to 300 minutes. The green contour is the detected contour and the black contour is the contour chosen to trace the light 
curves. 
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and then in the 193 channel. For the method 2 (phase 
difference), the time lag for 335-193 is inferior to the 
time lag for 335-171 (142 min compared to 155 min). 
The second method gives an ordering of channel that is 
more conform to the one predicted from a monotonic 
temperature decrease and from the AIA response func¬ 
tions. Only the 94-335 and 171-131 pairs give a time 
lag that is not expected in this predicted ordering. In¬ 
deed, the absolute value of the time lag of 94-335 (114 
min) is smaller than the one of 335-211 (122 min) and 
thus the intensity peaks in 94 before that in 211. Simi¬ 
lar ordering is found with method 1 as the time lags for 
94-335 (115 min) and 335-211 (113 min) are the same 
given the time resolution. For the 171-131 pair the 
time lag is -51 min with method 2 compared to -1 min 
with method 1. Method 2 (phase difference) is lim¬ 
ited in some passbands by the Fourier power spectrum 
at the detected frequency: as the normalized power is 
weaker for the 171 and the 131 channels, method 2 is 
more affected by noise and thus less reliable for the 
171-131 pair. 

3. Discussion 

We selected three cases of long-period intensity pul¬ 
sations and analyzed them in the six coronal chan¬ 
nels of AIA. These pulsations have periods of sev¬ 
eral hours, from 3.8 hours to 9.0 hours for the selected 
events. They are strongly visible in most of the pass- 
bands and are clearly matching loops structures in ac¬ 
tive regions. In each case, only some loops among 
those of the active regions are pulsating. 

3.1. Evolution of the thermal structure of pulsat¬ 
ing loops 

Multi-spectral analysis with AIA allows to proceed 
with a thermal study of these cases. The analysis of the 
thermal structure, i.e. the distribution and evolution 
of temperatures and densities, of the loops where pul¬ 
sations occur, reveals periodic variations. The DEM 
analysis also reveals high amplitude variations in the 
coolward wing of the DEM. We link these variations 
to a signature of heating of these loops. The variation 
of the DEM slope, which are large, can tell us about 
the timescales of the heating. Even though 
|et al.| ( |20T3| l have shown that the DEM slope 
constrained, we interpret here only the relative varia¬ 
tions of the slope. A steep slope implies that the heated 
plasma does not have enough time to cool down to 
lower temperatures, and then that the frequency of in¬ 


Guennou 
is poorly 


dividual heating events is high compared to the cool¬ 
ing time. On the contrary, when the slope is weak, a 
greater quantity of plasma has cooled down before be¬ 
ing re-heated, implying a lower frequency of heating. 

There are clear pulsations for the total emission 
measure for the three events studied. Eor event 1, 
there is a peak of normalized power above the adopted 
threshold (99% of confidence level) for the peak tem¬ 
perature. Eor event 2 and event 3, there are peaks of 
normalized power, respectively for the peak tempera¬ 
ture and the slope, below this threshold but still at high 
confidence levels. These small peaks have respectively 
97% (8.7cr) and 96% (8.4cr) of confidence level. Pul¬ 
sations do not come out clearly from the Eourier anal¬ 
ysis of the slope and the peak temperature. This can 
be explained by several effects impacting the DEM re¬ 
sults: 


Eirst, we remind that the DEM analysis is made 
without any subtraction of the background and fore¬ 
ground emissions because it would not be practical nor 
meaningful given that the detected pulsations are much 
longer than the typical life time of loops. However 
background emissions can be seen in images presented 
in Eig. Eig. 13 and Eig. 22 especially at 131 A, 
where we can clearly see the emission from the tran¬ 
sition region. Superimposition of structures with dif¬ 
ferent temperature and density conditions can partially 
hide the signal from the pulsating loops. 


Second, the importance of the background and fore¬ 
ground emission can be different in the different pass- 
bands. As an illustration, for event 1, where there 
are no pulsations (given the 99% confidence level) for 
the DEM slope, the normalized pulsating power is the 
strongest in the hotter channels (193, 211 and 335). 
The three other ones which could constrain the cool 
wing of the DEM may be the most affected by the 
background and foreground emission. On the contrary, 
for event 3, pulsations are exclusively visible in the 
cooler passbands: 94 (i.e. the cool component of 94), 
131, and 171. 


Third, Guennou et al. (2013 1 showed that the slope 
is poorly constrained and thus sensitive to random vari¬ 
ations in the light curves, which may explain why the 
time series of the slope is itself noisy. 


To conclude, the lack of clear pulsations (given the 
99% confidence level) for the DEM slope for the three 
events and for the DEM peak temperature for event 2 
and event 3, could come from a combination of the 
three above described effects. In addition even if it 
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does not come out clearly from the Fourier spectrum, 
there are indications of periodic pulsations in the DEM 
slope from correlations analysis. Therefore we con¬ 
clude that there are oscillations in the thermal structure 
of these loops. 


3.2. Widespread cooling 


The studied loops are seen in a state of cooling, as 
was shown using the same method as V&K12: time 
lag maps for six pairs of channels show widespread 
cooling for the three active regions, independently of 
the pulsating behavior of some of the loops in these ac¬ 
tive regions. To avoid the effect of the background and 
foreground emissions, we confirm this result using a 
second method, by extracting the pulsating component 
of the signal in Fourier space. 


But why do we only observe cooling? The usual 
explanation is that the density is too small for the 
plasma to be detectable during the heating phase. The 
widespread cooling can also be a signature of heating 
mechanisms: V&K12 and |Viall & Klimchuk] ( [20T3| ) 
interpret it as a signature of heating by nanoflare 
storms. In the nanoflare model indeed, loops strands 
are heated (by impulsive events) and cool indepen¬ 
dently (e.g Warren et al. 12003) . In that case, the cool¬ 
ing time of a bundle can be much longer than the 
typical life time of individual loops. Nevertheless, 
for Lionello et al. ( 2013| l, this widespread cooling is 
also consistent with quasi-continuous heating, as their 
model can reproduce long cooling times: they can be 
also interpreted by superimpositions of loops cooling 
independently along the line of sight and are not nec¬ 
essarily a signature of impulsive heating mechanisms. 


In the time lag maps we do not notice differences 
between the pulsating loops and the rest of the active 
regions. The cooling phase is thus the same for the 
entire active region but that does not exclude different 
processes of heating between the pulsating loops and 
the rest of the active regions. The heating mechanism 
may or may not be the same inside and outside of the 
pulsating region but there is not yet enough informa¬ 
tion to conclude. 


3.3. Possible physical explanations for long-period 
pulsations 

In order to explain long-period pulsations seen in the 
coronal EUV channels and in the thermal structure of 
warm loops, we can put forward at least two physical 
explanations depending of the supposed nature of the 


heating. If it is impulsive, we could envisage periodic 
nanoflare storms, i.e. nanoflare storms with periodical 
changes of the frequency of individual heating events, 
from high frequency (short time between nanoflare 
events) to a lower frequency (more time for the plasma 
to cool down before being re-heated). Then the DEM 
would change according to the changes in heating fre¬ 
quency. A change in the frequency of nanoflares might 
come from a coupling between the motions at the foot- 
points and the magnetic structure of the loops. 


On the other-hand, a quasi-continuous heating 
at the loops footpoints would also be a good can¬ 
didate to explain our observations. This scenario 
drives a condensation-evaporation cycle (thermal non¬ 
equilibrium): the plasma condenses at the top of the 
loops, this material falls toward the footpoints, then 
chromospheric material from the footpoints evapo¬ 
rates and fills the loops with hot plasma. Actually, 
even if the role of thermal non-equilibrium was of¬ 
ten ruled out in the heating of warm coronal loops 
(Klimchuk et al. 2010| l, we find that the long-period 
intensity pulsations observed in loops are consistent 
with condensation-evaporation cycles. Fig. [T^ repre¬ 
sents the time lag between the peak temperature and 
the emission measure for the entire field of view for 
event 1. This map is built by subtracting phase maps 
of the DEM peak temperature and of the emission 
measure at the main frequency detected for event 1. 
This allows us to obtain the time lag map between the 
peak temperature and the emission measure for only 
the periodic component of these parameters. There¬ 
fore, pulsating loops at the detected frequency are iso¬ 
lated in Fourier space. This map shows a time lag of 
123 min on average in the selected contour which is 
what we already found in Fig. without Fourier fil¬ 
tering. This corresponds to a de-phasing of 81°. This 


result is to be compared with Fig. 6 in Mikic et al. 


( 2013| l, where they compare the evolution of the tem¬ 
perature and the density at the apex of a symmetric 
simulated loop, with a nonuniform cross-section area 
and a nonuniform heating. The temperature and the 
density of their simulation have cyclic variations with 
a period of 5.0 hours, which is at the order of the pe¬ 
riods that we observed. The density curves have a 
delay of about 1.5 hours compared to the temperature 
curve. That implies that the temperature variations al¬ 
ways precede the density variations with a de-phasing 
of approximately 108°. For this simulation the authors 
note development of incomplete condensation cycles 
localized near the loop legs and a temperature at the 
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loop apex that remains above 1 MK. The observed and 


peak temperature-emission measure 


modeled de-phasings are comparable, as far as Mikic 
et al. ( 2013| l use the density and we have only access 
to the total emission measure. 


Even if long-period intensity pulsations are con¬ 
nected to heating processes in loops, these pulsations 
occur only in some loops of the studied active regions 
and not in the entire active regions. In other words, 
why don’t we have long-period intensity pulsations in 
all loops? This could be due to the background and 
foreground emission as developed earlier in the text. 
Intrinsic properties of loops could also explain why 
some loops do not show this kind of pulsations. For 
Mikic et al. ( 20I3| l, the plasma response to a steady 
heating mainly concentrated at the loops footpoints de¬ 
pends on the loops geometry: the uniformity and the 
symmetry of the cross-sectional area, the symmetry 
of the loop profile and the symmetry of the heating 
prohle with different heating functions. The geometry 
of the loops and of the heating could lead to different 
condensation strengths that they call ’’incomplete” and 
’’complete” condensation. This connects thermal non¬ 
equilibrium processes to both coronal rain and heating 
of warm loops. 

The wide distribution of observed periods could be 
due to inherent properties of the loops. In lAuchere 


et al.| ( ^I4| l a possible relation between the length 
of the loops and the period of the pulsations was dis¬ 
cussed. Here also we notice that for these three events, 
the longer the period of pulsations, the larger the de¬ 
tection area. It remains to be confirmed that the loop 
length is a good candidate to explain the different pe¬ 
riods observed in pulsating loops. 

4. Summary and Conclusions 


Auchere et al. P014|l have shown that long-period 


intensity pulsations are a common phenomenon in 
coronal loops. Using the same detection algorithm, we 
have detected many other cases in the six coronal chan¬ 
nels of AIA. Among these, we selected three events de¬ 
tected in three different active regions. The observed 
pulsations have respectively periods of 9.0 hours, 5.6 
hours, and 3.8 hours. The normalized power maps of 
the three studied active regions clearly show that these 
pulsations are strongly localized in some loops, but not 
in all loops, of active regions. Thus, either the occur¬ 
rence of these cycles depends on specific loop prop¬ 
erties or some events could be masked by the back- 
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Fig. 12.— Time lag map obtained from differences of 
phase maps of the peak temperature and of the emis¬ 
sion measure. We explored time shifts from -300 min¬ 
utes to 300 minutes. The green contour is the detected 
contour and the black contour is the contour manually 
selected. 

ground and foreground emissions. 

We investigate the physical properties underlying 
these pulsations. First by means of a DEM analy¬ 
sis using the active region DEM model from |Guen-| 
nou et al. (20131. We notice periodic variations of 


the shape of the DEM during the sequences studied. 
The thermal structure of these loops changes period¬ 
ically from a strong multithermal structure to a more 
isothermal structure. We link these variations to a 
signature of heating mechanisms. Moreover, using 
the same method than in V&K12, we notice that the 
studied active regions are seen mainly in their cool¬ 
ing phase. This is confirmed by extracting in Fourier 
space the loops at the detected frequency, allowing to 
remove contributions from the background and fore¬ 
ground emissions. Although this widespread cooling 
could be linked to nanoflare heating, it does not in¬ 
evitably imply impulsive heating. In fact, Lionello 


et al. ( 2013| l shows, by means of 3D-hydrodynamic 
simulations, that thermal non-equilibrium processes 
could be also consistent with widespread cooling. Our 
observations bring a new element to this debate. Since 
these pulsations appear to be common in coronal loops 
and are linked, as we suggest, to the heating of the 
plasma, then any model of loop heating must be able 
to reproduce these new observations. 
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The SDO/AIA images are available by courtesy of 
NASA/SDO and the AIA science teams. This work 
made extensive use of the AIA archive at MEDOC, 
http://medoc-sdo.ias.u-psud.fr 
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A. Event 2 
A.I. Data sample 


Fig. [T^ shows images of a second exemple of pulsating loops with long-period on December 30, 2012 at 21:28 UT. 
For this event 2, the NOAA AR 11637 was tracked from December 28, 2012 10:00 UT to January 02, 2013 18:42 UT 
(i.e. about 129 hours). As for event 1, we focused on the three middle days of this sequence, from December 29, 
2012 14:26 UT to January 01, 2013 14:29 UT (i.e. 72 hours). The contour in green is the contour detected in the 
171 A channel centered at 342.5° of longitude and 20.5° of latitude with an area of 12.9 heliographic square degrees 
(1495 Mm^). We manually selected a small contour (in black) included in the detected contour, to obtain the time 


series presented in Fig. 14 


We display the normalized Fourier spectra in Fig. 15 In the 94 A, 131 A, 171 A, 193 A, and 335 A channels there 
is a significant peak of power at 49.9 juHz (5.6 hours). The normalized powers are respectively 16.7cr, 13.5cr, 19.9cr, 
19.3cr, and 25.5cr. Even if there is no peak above lOcr in 211, we can notice a small peak with a normalized power 
of 9.6cr at 49.9 juHz. At 94 A, and 193 A there is also a small peak at 42.2 //Hz with respectively 8.1cr and 11.7cr of 
normalized power. 


Fig. [T^represents the maps of normalized power for event 2, at the detection frequency (49.9 //Hz, i.e. 5.6 hours). 
Loops of NOAA AR 11637 are clearly seen in these maps for all the passbands. The averaged normalized power is 
larger than lOcr for all the passbands, except 211 A (7cr). Loops of NOAA AR 11640, a second active region around 
320° in longitude and 30° in latitude, have also normalized power above lOcr. The automatic algorithm does not detect 
these loops with the set of arbitrary thresholds. This suggests that more of these events are probably not detected by 
our code. Loops of this second active region (emerging at the beginning of the sequence) seems to have almost the 
same length, that can be a reason why they have the same frequency of pulsations (see section [33]l. 


A.2. Differential emission measure (DEM) analysis 


The evolution of the mean DEM slope, peak temperature, emission measure, and in the selected contour of 
event 2 is presented in Lig. 14 Lig. [TT] represents the power spectra of these parameters. Lor the emission measure 
there is a peak power at 49.9 //Hz, i.e. 5.6 hours. The power at the peak is 15.0cr. If we look under the detection 
threshold at 49.9 //Hz, there is a peak with a power of 8.0cr for^^, but its variations from to 0.3 to 1.8 indicates a good 
fit. Lor the peak temperature, there is a small peak (8.7cr) at 46.1 //Hz. Lor the slope there is no significant peak of 
power. 


Even if the power is weak for the peak temperature and there are no pulsations found with the Lourier analysis for 
the slope, we found a clear correlation between the slope and both the peak temperature and the emission measure, 
between 97 hours and 120 hours after the beginning of the full sequence. In Lig.[^ we plot the cross-correlation values 
for DEM peak temperature vs. emission measure and slope. We find a clear correlation between the peak temperature 
and the emission measure (cross-correlation value of 0.66) and anti-correlation between the peak temperature and the 
slope (cross-correlation value of -0.74). We can conclude that there are indications that the slope varies periodically 
with periods of hours and thus that the thermal structure of the studied loops undergoes periodical changes. 


The amplitudes of variations of the DEM parameters are larger than for event 1: there are variations from 2.2 to 
3.2 MK (a relative amplitude of 45%) for peak temperature, from 3.2 x 10^^ to 6.4 x 10^^ cm“^ for emission measure 
(a relative amplitude of 97%), and from 1.7 to 4.8 for the slope (a relative amplitude of 176%). 


A.3. Evidence for widespread cooling 



pairs , mainly negative for 94-335 and equal to zero for 171-131. The mean time lags in the black contour can be 
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Fig. 13.— Same as Fig.[^for event 2 on December 30, 2012 at 21:28 UT. This event is localized in NOAA AR 11637. 
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Fourier power spectrum (94 A) 
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Fig. 15.— Fourier power spectra for event 2, same as Fig.[^ 


23 





















































Latitude (degrees) Latitude (degrees) 



320 330 340 350 

Longitude (degrees) 


S’ 

5 , 


■D 

D 

3 





193 A 


320 330 340 350 

Longitude (degrees) 


320 330 340 350 

Longitude (degrees) 


335 A 


320 330 340 350 

Longitude (degrees) 



320 330 340 350 

Longitude (degrees) 


320 330 340 350 

Longitude (degrees) 



Fig. 16.— Same as Fig.|^with the same scale. Maps of normalized power at the frequency of detection (49.9 pHz, i.e. 
5.6 hours) for event 2. The black contour, included in the detected contour at 171 A (in green), is the contour manually 
selected for detailed time series. 
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Fig. 17.— Fourier power spectra for the three parameters of the DEM model for event 2, same as Fig.|^ 



Time shift (minutes) 


Fig. 18.— Same as Fig. We used the averaged curves at 1 min of cadence averaged over the small black contour 
presented in Fig.[2between January 03, 2013 22:28 UT (97 h after the beginning of the long sequence) and January 
04, 2013 21:28 UT (120 h after the beginning of the long sequence). 
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Fig. 19.— Same as Fig.j^for event 2 with the same scale. 
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compared in Table Using the cross-correlation method, intensity peaks in the following channel order: 335, 211, 
171 and 131, 94, and 193. Using the phase difference method, the order becomes: 335, 211, 193, 94, 171, and 131, 
that is the order expected from the expected temperature evolution and from the AIA response functions. This second 
method is better in this case (event 2) because the power is high enough in all the passbands at the detected frequency. 
However, for pairs like 94-335 and 335-193 time lags are almost equal. We can notice that these time lags are smaller 
than the time lags found for event 1 (about the half). These smaller time lags are likely to be due to the smaller loop 
length (smaller detection area) for event 2. 


B. Event 3 
B.l. Data sample 


Fig. 22 shows images of a last example, event 3, on August 09, 2011 12:12 UT. For this event, the ROI (45° in 
longitude) was tracked from August 07, 2011 04:00 UT to August 13, 2011 05:49 UT (i.e. about 146 hours), and we 
select the three middle days of this sequence, from August 08, 2011 16:53 UT to August 11, 2011 16:56 UT (i.e. 72 
hours). Pulsations are detected in loops located above the core of NOAA AR 11268. The detected contour at 94 A (in 
green) is centered on 218.0° of longitude and 21.4° of latitude with an area of 6.3 heliographic square degrees (730 
Mm^). In FigU^we represent the average light curves in the selected contour (in black). The corvesponding Fourier 
spectra in Fig^j^show strong signal in the 94 A, 131 A, and 171 A channels. Pulsations are detected in these three 
bands at 72.9 /rHz (3.8 hours) with respectively a power of 27.4cr, 18.4cr, and 14.7cr above the estimated background 
spectrum. At 211 A and 335 A there is no significant power peak, but at 193 A there is a small peak at 8.2cr. In this 
last channel, there is also a small peak (7.0cr) at 88.3 yuHz. 


Fig. [^represents the maps of normalized power for event 3, at the detection frequency (72.9 /iHz, i.e. 3.8 hours). 
The normalized power in the 94 A, 131 A, and 171 A channels is higher than 14cr on average inside the black contour. 
For the other passbands, the normalized power is weaker, but the power maps also display loops near the core of the 
active region. 


B.2. Differential emission measure (DEM) analysis 


The evolution of the mean DEM slope, peak temperature, emission measure, and in the selected contour for 


event 3 is presented in Fig. 23 


Fig, [^represents the power spectra of these parameters. Pulsations are detected for the emission mesure and;^f^ 
at 72.9 yuHz, the same frequency of detected pulsations in the light curves. The normalized power at these peaks are 
respectively 11 .Ocr and 10. Icr. If we look under the detection threshold, there is a small peak of power for the the slope 
(8.4cr). Between 40 hours and 70 hours after the beginning of the large sequence we found a clear correlation (0.78 of 


peak cross-correlation value) between the peak temperature and the emission measure as seen in Fig. 27 


The amplitudes of variations of the DEM parameters are larger than for event 1 and event 2: there are variations 
from 1.9 to 3.8 MK (a relative amplitude of 99%) for the peak temperature, from 3.0 x 10^^ to 8.9 x 10^^ cm“^ for 
the emission measure (a relative amplitude of 196%), and from 1.4 to 3.3 for the slope (a relative amplitude of 144%). 


Table 2: Sa me as TablefTIfor event 2. 


Pairs of channels 


335-211 

211-193 

335-193 

94-335 

335-171 

171-131 


Time lag (min) 
from cross-correlations (Pig. 


19i 


35 

25 

80 


-71 


69 

1 


Time lag (min) 


from differences of phase Pig .(21 1 


37 

31 

68 


-68 


73 

-4 
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Fig. 20.— Same as Fig. 10 for event 2 with the same scale. 
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Fig. 21.— Same as Fig. 11 for event 2 with the same scale. 


29 












Latitude (degrees) Latitude (degrees) 



205 210 215 220 225 230 

Longitude (degrees) 


131 A 171 A 



205 210 215 220 225 230 205 210 215 220 225 230 

Longitude (degrees) Longitude (degrees) 


193 A 211 A 335 A 



205 210 215 220 225 230 

Longitude (degrees) 


205 210 215 220 225 230 

Longitude (degrees) 


205 210 215 220 225 230 

Longitude (degrees) 


Fig. 22.— Same as Fig. [T]for event 3 on August 09, 2011 12:12 UT. This event is localized in NOAA AR 11268. 
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Fig. 24.— Fourier power spectra for event 3, same as Fig.[^ 
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Fig. 25.— Same as Fig.|^with the same scale. Maps of normalized power at the frequency of detection (72.9 pHz, i.e. 
3.8 hours) for event 3. The black contour, included in the detected contour at 94 A (in green), is the contour manually 
selected for detailed time series. 
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These strong variations of both the peak temperature and the total emission measure are observed in the core of an 
active region where any type of activity is in general expected to be the highest. 


B.3. Evidence for widespread cooling 

The time lags maps in Fig. (made with the peak cross-correlation method) highlight a widespread cooling in this 
active region. Fig. [^represents the phase maps. The phases are constant along the loops. Maps produced with the 
second method (phase differences), shown in Fig. 30 lead to the same conclusion of widespread cooling. The delay 
presented in Table for both methods, show a better succession of channels with the second method; with the first 
method we have the 335, 211, 94, 193, 171, and 131 order, whereas with the second we obtain 335, 211,193 and 94, 
171 and 131, which is the ordering expected for the AIA channels. However for this event and for both methods, 
results are noisy for some pairs of channels. Indeed, for the first method, the peak cross-correlation values are poor 
for some pairs in the black contour (Fig. |28l l and the time lags for 335-193 and 335-171 are the most affected. The 
second method suffers from the lack of strong power in the 335, 211 and 193 channels. For event 3, we observe time 
lags smaller than for event 1 and event 2, that can be explained by the probably smaller loop length (smaller detection 
area) for event 3. 


Table 3; Sa me as Table [T] for event 3. 


Pairs of channels 


Time lag (min) 
from cross-correlations (Fig. 


Time lag (min) 


28 I from differences of phase Fig .(30 1 


335-211 

26 

4 

211-193 

14 

29 

335-193 

50 

33 

94-335 

-36 

-33 

335-171 

76 

42 

171-131 

0 

-1 
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Fourier power spectrum (peak temperature) 



Fourier power spectrum (emission measure) 



Fourier power spectrum (siope) 



Fourier power spectrum (x^) 



Fig. 26.— Fourier power spectra for the three parameters of the DEM model for event 3, same as Fig.|^ 



Time shift (minutes) 


Fig. 27.— Same as Fig. We used the averaged curves at 1 min of cadence averaged over the small black contour 
presented in Fig.[T]between August 09, 2011 00:00 UT (40 h after the beginning of the long sequence) and August 10, 
2011 06:00 UT (70 h after the beginning of the long sequence). 
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Fig. 28.— Same as Fig.j^for event 3 with the same scale. 
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Fig. 29.— Same as Fig. 10 for event 3 with the same scale. 
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Fig. 30.— Same as Fig. 11 for event 3 with the same scale. 
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